From modes to movement in the behavior of C. elegans 
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Organisms move through the world by changing their shape, and here we explore the mapping 
from shape space to movements in the nematode C. elegans as it crawls on a planar agar surface. 
We characterize the statistics of the trajectories through the correlation functions of the orientation 
angular velocity, orientation angle and the mean-squared displacement, and we find that the loss 
of orientational memory has significant contributions from both abrupt, large amplitude turning 
events and the continuous dynamics between these events. Further, we demonstrate long-time 
persistence of orientational memory in the intervals between abrupt turns. Building on recent 
work demonstrating that C. elegans movements are restricted to a low-dimensional shape space, we 
construct a map from the dynamics in this shape space to the trajectory of the worm along the 
agar. We use this connection to illustrate that changes in the continuous dynamics reveal subtle 
differences in movement strategy that occur among mutants defective in two classes of dopamine 
receptors. 



I. INTRODUCTION 

From the swimming motions of E. coli pQ to the mobil- 
ity of human populations [2] , the way in which organisms 
move through the world profoundly influences their expe- 
rience. Ultimately, these strategies of movement change 
the chances for survival and reproduction, and thus are 
subject to natural selection. Crucially, whether a move- 
ment is adaptive depends on its relation to the outside 
world, but the organism has control only over its inter- 
nal states. From this internal point of view, movement is 
not translation or rotation of the body relative to a fixed 
external coordinate system, but rather transformations 
of the organism's shape as measured in its own intrinsic 
coordinates. 

In general, the connection between transformations in 
shape space and movement through the world is compli- 
cated. There is a long tradition of work which tries to 
make this connection through analytic approximations of 
the equations describing the mechanics of the organism's 
interaction with the outside world. This approach is per- 
haps best developed for swimming and flying organisms 
[3] [4], and there are particularly elegant results in the 
limit of swimming at low Reynolds' number [5j|6j[7]. All 
of these methods depend on some small parameter in the 
physical interaction between the organism and its envi- 
ronment. A very different possibility for simplifying the 
relation between shapes and movement arises if the space 
of shapes itself is limited. In several systems, the poten- 
tially high dimensional space of shapes or movements is 
not sampled uniformly under natural conditions, so that 
one can recognize a lower dimensional manifold that fully 
describes the system [8j [9] [TUJ [11]. In these cases it is 
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possible to ask empirically how motions on this low di- 
mensional manifold map into movements relative to the 
outside world. 

The motion of the nematode C. elegans provides an 
example of this dimensionality reduction. In previous 
work, we found that the shapes taken on by the worm's 
body are well-approximated by a four dimensional space 
spanned by elementary shapes or 'eigenworms' [12] . Here 
we connect the dynamics in this low dimensional space 
of shapes to the trajectories of worms as they crawl on 
an agar plate, the conventional experimental setup for 
studying worm behavior [13] [14] [15] [16]. In the process, 
we offer a new analysis of the trajectories themselves, 
and show how the intrinsic shape dynamics gives us a 
more refined tool for the analysis of mutant locomotory 
behavior. 



II. CENTER OF MASS TRAJECTORIES 

To understand how the motion of the worm is deter- 
mined by its internal dynamics in shape space, we first 
characterize the motion itself. In Fig [lji, we show an 
example of the worm's trajectory, defined as the cen- 
troid of the black and white body image. These data 
are collected as in Ref [12], using tracking microscopy, 
and what is shown is one example from a data set that 
includes TV = 33 wild- type worms, each tracked for 35 
minutes. The worms were transfered to the agar plate 
using a platinum worm pick and we excluded the first 
400 seconds of each tracking run to avoid any influences 
from this mechanical stimulation. 

The trajectory consists of gently curving segments, in- 
terrupted by sharp turns, often associated with charac- 
teristically curved body shapes; these discrete events are 
known as Sl-turns. It is therefore tempting to think of 
these trajectories as being approximately like those of E. 
coli, consisting of long, relatively straight runs punctu- 
ated by tumbles, which randomly reorient the cell [T7] . 
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FIG. 1: Centroid trajectories, (a) A typical center of mass 
trajectory on the agar plate. The path includes both abrupt 
orientation changes (red) in which the worm shape is deeply 
bent (inset), and continuously curving segments, (b) Defini- 
tion of the local orientation angle a(t) and the local orienta- 
tion velocity u(t) = da(t)/dt. The inset shows the cumulative 
distribution of orientation velocities (red) contrasted with a 
Gaussian distribution that has the same variance (blue) . The 
centroid time series {x cm (t), y C m(t)} was filtered with a third- 
order polynomial in a running window 11 frames in length 
and a cm and co cm were built using this filter. 



Indeed, variations of this model in which trajectories are 
segmented into discrete 'runs' have long been used in 
study of C. elegans [T3j HE] and in the analysis of move- 
ments for a wide variety of organisms p~9j [20] . 

Trajectories are characterized by the speed of motion, 
v C m = \x\, the local tangent angle, a(i), and the rate 
of local curvature uo = da/dt at each moment in time 
(Fig [IJ)). The standard deviation of the curvature is 
fo rms = 0.223 ±0.001 rad/sec, but p(uo) has long tails, as 
shown by comparing the cumulative distribution P(\oo\) 
to a Gaussian distribution with the same variance (inset 
to Fig[^>). Excursions to the large amplitude tails corre- 
spond to abrupt reorientation events and are colored red 
in Fig la. 

If worms lose orientational memory then their move- 
ments will look diffusive at long times. One signature of 
such behavior is contained in the mean-square distance 
between two points on the trajectory, 



(( Sx ) 2 ) T = (\x(t+T)-m 2 ), 



(i) 



which for diffusion should grow linearly with the time r. 
In Fig [2^i, we see that this is what happens for times 
longer than r ~ 10 s. At shorter times, the mean-square 
displacement grows as the square of the time difference, 
which corresponds to ballistic motion at a fixed velocity. 
This suggests that the time scale for worms to lose di- 
rectional memory is ~ 10 s, and this can be seen directly 
from the correlation function for the local tangent angle, 



C a (r) = (cos[a(t + r) -a(t)]), 



(2) 



as shown in Fig(2}}. 

A sequence of uncorrelated runs and tumbles trans- 
lates into a mathematical model in which the orientation 
angle a is, from time to time, completely randomized by 
discrete events. In E. coli these events are the tumbles, 



and it is natural to think that in C. elegans these events 
are the Q-turns. If turns occur randomly at rate r tU rn, 
and generate completely new, uncorrelated directions of 
movement, then the correlation function for direction will 
decay as 



C a (r) = exp(-r tU rn|r| 



(3) 



But since r tU rn typically is less than two per minute 
[T3l H8] . Q-turns alone can't explain the shorter time 
of decay of the directional correlations, as seen quantita- 
tively in Fig[l]i. Thus, changes in the orientation angle 
a occurring in between ^-turns must play a key role. 

To find an objective but tractable definition of the dis- 
crete events that interrupt the trajectory, we examine the 
approximations made in drawing the trajectory itself, as 
in Fig [I] Following general practice, we generated tra- 
jectories by tracking the center of mass of the worm's 
silhouette, x cm (t). An alternative would be to follow the 
midpoint of the worm as indicated in Fig|3^i, generating a 
trajectory x m id(t). When the worm moves smoothly and 
doesn't change shape dramatically, these different defini- 
tions of trajectory produce nearly identical results. But, 
during those moments when the worm pauses and reori- 
ents, these trajectories diverge. This gives us a natural 
way to identify abrupt reorientation events. 

We look at the difference between the center of mass 
and the midpoint trajectories, Sx = x cm — £ m id and de- 
fine the difference speed v$ = Sx. From the conditional 
density p(vs\v cm ), shown in Fig [3}}, we can see that this 
difference speed v§ is typically much smaller than the 
center of mass speed v cm . But at small v cm there is sub- 
stantial probability that v§ > v cm and when this occurs 
the centroid orientation is no longer an accurate measure 
of the worm's heading. At low speeds there is also a large 
variance in local curvature, which can be seen in the con- 
ditional density p(^ C m|^cm), Fig [3^. Thus our segmenta- 
tion based on v cm is similar to those based on large values 
of co cm . To simultaneously isolate the largest amplitude 
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FIG. 2: Trajectory correlations, (a) Mean-square displace- 
ment as a function of time along the trajectory [Eq Q]. At 
short times we see ballistic (Sx oc r) behavior, crossing over 
to diffusion [(Sx) 2 oc r] at long times, (b) The orientation an- 
gle correlation function [Eq |2])] for worm motion (red), con- 
trasted with a simple "run and tumble" model (blue) in which 
reorientation events occur randomly with rate m — 0.045 Hz, 
corresponding to the rate of Q-turns [131 [T8] . 
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turns and maintain the validity of the centroid orien- 
tation within intervals we choose ^thresh = 0.035 mm/s. 
Since lower centroid speeds are associated with curvier 
shapes, the choice of a threshold in v cm amounts to tol- 
erating a fraction of frames with v§ > thresh and this 
fraction is 25% at our threshold. Most of the time v cm is 
larger and the fraction substantially smaller. 
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FIG. 4: The loss of orientation memory is captured by inde- 
pendent contributions from abrupt reorientation events and 
continuous dynamics between these events. The orientation 
correlation function CJ, nterval (r) computed by averaging over 
all delays within the same interval is shown in green. The 
probability P(r) that two times separated by r are within 
the same interval is shown in blue. The predicted orientation 
correlation function from Eq Q is shown in black, compared 
with the data, shown as red points with standard errors of 
the mean. 



FIG. 3: The reorientation of directional heading, (a) The 
displacement vector Sx = x cm — x m id- (b) The conditional 
probability density p(v$\v cm )- At low speeds, the centroid ve- 
locity mixes changes of the body shape with changes in the 
directional heading, (c) Low centroid speeds are also asso- 
ciated with large centroid orientation changes as shown by 
the conditional distribution p((jj\v cm ). The solid white line 
denotes the speed threshold at which we separate continuous 
intervals from discrete events. 

Using the intervals between reorientation events as de- 
fined above, we can disentangle the contributions of con- 
tinuous and discrete reorientation. We compute the ori- 
entation correlation function C a (r) within the intervals, 
and the result, C^ nterval (r), is shown as the green curve 
of Fig. [4] Removing the intervening abrupt reorienta- 
tion events reveals angular correlations which decay on 
a much longer time scale. Furthermore we observe anti- 
correlations which indicate a marked departure from or- 
dinary orientational diffusion. 

As an approximation, we assume that the continu- 
ous dynamics of orientation within arcing segments of 
the trajectory are independent of the discrete reorienta- 
tion events, and that the reorientation events are instata- 
neous. Then the full orientation correlation function is 
the product of two terms, 



/ototal/^N r 1 



interval 



(r) x P( T ), 



(4) 



where P(r) is the measured probability that two points 
separated by time r come from the same interval. Figure 
[4] shows that this is a excellent approximation, with no 
adjustable parameters. 



III. MODES AND MOTIONS 

Following the center of mass of an organism is a con- 
venient but approximate measure of its movement strat- 
egy. A fuller, and organism centered, representation of 
locomotion is contained in the shape of the worm itself. 
In previous work [12 we found that the space of shapes 
of C. elegans during free locomotion is low-dimensional 
with four principal dimensions (eigenworms) capturing 
approximately 95% of the variance of the space of shapes. 
A summary of these results is shown in Fig [5] 

To construct the connection between modes and move- 
ment we confine our analysis to continuous intervals, 
within which co crn is a faithful measure of the worm's 
orientational velocity. We then try to construct a map 
from the amplitudes of motion along the four intrinsic 
shape modes into uo(t). The simplest model is a linear 
one, 



^modes 



(*) 



i=l 



(5) 



and by adjusting the coefficients f3[ we can capture 55% of 
the variance in uo. More generally we can try a nonlinear 
relationship, and allow that the current curvature of the 
trajectory also depends on the time derivatives of the 
shape. To do this, we introduce twelve variables, 



Zi=l:4 = 01:4 CO, 

da\.^{t) 



^i=5:J 



dt 



(6) 
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FIG. 5: Low dimensional description of worm shape, follow- 
ing Ref [T2] . (a) Each raw image of the worm is processed by 
passing a curve through the center of the body; the red circle 
marks the head. Arc length s along this curve is normalized, 
< s < 1, and we define the tangent t(s) = ^| to the cen- 
terline curve, (b) We rotate all images so that (0) is zero and 
thus 0(s) provides a description of the worms shape that is in- 
trinsic to the worm itself, (c) We decompose each shape 0(s) 
into contributions from the leading four eigenmodes, which 
capture 95% of the variances in shape space, and the am- 
plitudes of fluctuation along each mode are normalized to 
unit variance, (a?) = 1. (d) Fluctuations along the first two 
modes correspond to an oscillation, or nearly circular orbit 
in this projection of the shape space, (e) Fluctuations along 
the third mode show strongly non-Gaussian statistics. The 
first two modes drive the propulsive wave along the body, the 
third mode makes the dominant contribution to turning, and 
the fourth mode (not shown) describes localized fluctuations 
of the head and tail. 



d 2 ai. A (t) 
dt 2 ' 



and construct an estimate of the curvature as 

^modes 



ijk 



(8) 



(9) 



1J 



As shown in Fig |6^l, this approach captures more than 
80% of the variance in the curvature. More importantly, 
we also capture the dynamics of the curvature, as shown 
through the correlation function 



i 



(uj 2 



■((j(t + T)Lj(t)} 



(10) 



in Fig [6)3. 

The curvature is dominated by contributions from the 
amplitude of the third mode, which plausibly is con- 
nected to bending of the worm's body [12] . We have 
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FIG. 6: The intrinsic dynamics of worm shape capture the 
extrinsic motions of foraging trajectories, (a) We construct 
an eigenworm model of the orientation velocity (Eq. [9| and 
show the joint distribution p(cj cm , w mo des)- (b) The orienta- 
tion velocity two-point correlation function for the eigenworm 
model (blue) and the centroid motion (red). 



shown previously that the instantaneous speed of the 
worm v cm (t) oc 0(t), where the phase <p as the angle in the 
plane {ai,<22} (cf Fig[5]i). Taken together, then, we can 
map from the modes to the speed and local curvature 
of the center of mass trajectory, so we can completely 
reconstruct the worm's movements from its shape as a 
function of time. 



IV. MUTANTS AND ADAPTATION 

As with all organisms, C. elegans behavior is mod- 
ulated by the organism's experience. As an example, 
the rate of ^-turns decreases systematically with time 
away from food [13], as well as changing in response to 
thermal [21] and chemical [18] stimuli. It is also known 
that dopamine plays a significant role in these behaviors 
[22] [23] . Here we use the analytic tools developed above 
to characterize both adaptation and the dopamine mu- 
tants dop-2 (vsl05) and dop-3 (vsl06). 

To make connections with earlier work, we follow 
Ref [18] and identify Q-turns as moments when co cm > 
0.87 rad/s. With this threshold, we can measure the cu- 
mulative distribution of inter-turn intervals, or equiva- 
lently the probability that a turn has not yet occurred 
after a time r, and this is shown in Fig [7] During the 
course of our experiments the worms spend 40 min away 
from food; allowing for an initial adjustment to being 
placed on the the plate we divide the last 28.3 min into 
three equal epochs to search for adaptation to the envi- 
ronment. We see that the interval distributions vary sys- 
tematically with time. As reported previously, the mean 
time between turns gets longer as the worms spend more 
time away from food. In more detail, we see that the 
distribution has two components, 

P(r) = ae~ rsWtr + (1 - a)e- rio ^ r , (11) 

and only the slower component ri ong contributes to the 
lengthening of the times between turns. Repeating the 
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analysis for the dopamine mutants, we again find that the 
short time behavior of the distribution, summarized by 
^short is unchanging, while the long time behavior varies 
across time (Table 1). We also note that differences in 
Sl-turns between the dopamine mutants are statistically 
insignificant though both mutants turn more frequently 
than wild- type worms. 




FIG. 7: The cumulative distribution of times between Q-turns 
in wild-type worms. Q-events are denned as described in the 
text, and dashed lines denote bootstrap standard deviations. 



also combine though Eq ([9| to produce different turning 
dynamics and, in particular, dop-3 animals make longer- 
lasting turns which result in curvier trajectories. While 
these pheonotypic differences are subtle and have not 
been reported before, they are immediately apparent in 
the dynamics in shape space. Taken together we find 
continuous motion, including gradual re-orientation, and 
discrete turning behaviors are under independent genetic 
control. 
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FIG. 8: Correlation functions for the shape modes a\{t) and 
a3 (t), comparing the wild type N2 with the two mutant 
strains dop-2 and dop-3. Shorter timescales for the dop-3 mu- 
tants results in faster centroid speed and curvier trajectories. 
Dashed lines denote bootstrap standard deviations. 





r short (Hz, ±0.02) 


n ong (Hz, ±0.002) 


genotype 


early 


middle 


late 


early 


middle 


late 


N2 


0.63 


0.76 


0.63 


0.053 


0.040 


0.027 


dop-2 


0.68 


0.67 


0.68 


0.072 


0.048 


0.044 


dop-3 


0.65 


0.66 


0.60 


0.070 


0.046 


0.042 



TABLE I: Properties of the inter-turn interval distribution. 
Data as in Fig [7] were collected also for N = 50 dopamine mu- 
tants, eac h al so observed for 35 minutes, and all the data were 
fit to Eq |TT| ). Results are shown for the two rates, r s hort and 
Hong, that define the time scales for turning. We note that 
Rshort = 0.66 ± 0.04 Hz across all epochs and genetic variants, 
and that differences in Q-turns between the dopamine mu- 
tants are statistically insignificant. 

Our analysis above shows that continuous re- 
orientation in between turns is a significant component 
of the worm's motion, and that this behavior is driven 
by the dynamics in shape space. Indeed, although the 
statistics of turning are the same for the two different 
mutants, Fig [8] shows that the dynamics along the differ- 
ent modes are in fact quite different. Along modes 1 and 
2 (which form a quadrature pair), dop-2 is similar to the 
wild type, but dop-3 exhibits a faster oscillation. Recall- 
ing that this corresponds to the undulatory wave along 
the body [12], we predict that the dop-3 mutant should 
move more quickly, and this is observed when we track 
the center of mass motions: mean speeds of the three 
variants are v cm = [0.076,0.082,0.121] ± 0.001 mm/s for 
N2, dop-2 and dop-3 respectively. The mode dynamics 



V. DISCUSSION 

It is tempting to think of behavior as a sequence of 
discrete actions, each taken in consequence of a specific 
decision. In the simplest cases, such as the running and 
tumbling of E. coli, we can see these discrete events re- 
flected directly in the trajectory of the organism's cen- 
ter of mass motion pQ. In contrast, we find that, for 
C. elegans, discrete behaviors are roughly only half of 
the story: the exploratory trajectories of the worm get 
approximately equal contributions from discrete turning 
events and from continuous re-orientational motions in 
between the turns. Further, we can trace these contin- 
uous motions back to the underlying dynamics in the 
space of body shapes. Finally, we see that these different 
components of the motion are under independent genetic 
and adaptive control. 

Quantitatively, we found that the exploratory motions 
of C. elegans are composed of two principle reorienta- 
tion elements: abrupt events, including classical Sl-turns, 
which occur at infrequently, and the continuous dynamics 
of orientation between these events. These two processes 
both make significant contributions to the worm's total 
loss of orientational memory on the ~ 10 s time scale. 
By focusing on the continuous intervals between discrete 
turns, however, we find that the worm's orientation can 
exhibit a longer term memory, lasting two minutes or 
more, and anti-correlations, corresponding to an abun- 
dance of arcs as seen, for example, in the trajectory of 
Fig. la. The presence of these arcs as well as the differ- 
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ential behavior of the mutants suggests that C. elegans 
controls more aspects of its motion than the stochastic 
rate of abrupt events. Evidence for this sort of contin- 
uous control has also recently been observed during C. 
elegans chemotaxis [24] . 

Extending previous work showing that the body shapes 
of C. elegans are captured by four modes [12 j, we built 
a phenomenological model that connects the intrinsic 
dynamics of these modes to the speed and curvature 
of the worm's trajectory through the external world. 
This model allows us to connect the body configurations, 
which are what the neuromuscular system can control, to 
the behaviors that have adaptive value. As an example 
of what can be learned from this analysis, we studied the 
motion of two mutants, dop-2 and dop-3, which contain 
defective receptors for the neurotransmitter dopamine, 
an important component in the modulation of foraging 



strategy. Although these mutants have nearly identical 
statistics when we look at their discrete turning behav- 
iors, their continuous motions, as seen in the dynamics of 
fluctuations along four different shape dimensions, show 
substantial differences. This suggests that the goal of 
mapping genes to behavior [25j [26] will require us to look 
much more closely, and quantitatively, at the behavior of 
individual organisms. 
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